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AERONAUTIC SYMBOLS 
I. FUNDAMENTAL AND DERIVED UNITS 




w 

Q 

m 

I 

v 

S 

S„ 

o 

b 

c 

A 

V 

i z 
L 
D 
Do 


D t 

D P 



Symbol 

Metric 

English 

Unit 

Abbrevia- 

tion 

Unit 

Abbrevia- 

tion 

Length 

Time 

Force 

l 

t 

F 

meter 

second 

weight of 1 kilogram 

m 

s 

kg 

foot (or mile) _ 

second (or hour) 

weight of 1 pound 

ft (or mi) 
sec (or hr) 
lb 

PoWPT 

P 

V 

horsepower (metric) 


horsepower 

hp , 

mph 

fps 

I UW vl w — — — 

Speed 

/kilometers per hour 

(meters per second 

kph 

mps 

miles per hour 

feet per second 


2. GENERAL SYMBOLS 


Weight = mg 

Standard acceleration of gravity =9. 80665 m/s 3 
or 32.1740 ft/sec 2 
W 

Mass=— 

g 

Moment of inertia = mk*. (Indicate axis of 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


v Kinematic viscosity 

p Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m'^-s 2 at 15° C 
and 760 mm; or 0.002378 lb-ft“ 4 sec 2 
Specific weight of “standard” air, 1,2255 kg/m 3 or 
0.07651 lb/cu ft 


3. AERODYNAMIC SYMBOLS 


Area 

Area of wing 
Gap 
Span 
Chord 




Aspect ratio, g 
True air speed 
D 3 mamic pressure, ^pT 73 

Lift, absolute coefficient C L =^ 


D 




Drag, absolute coefficient 
Profile drag, absolute coefficient 
Induced drag, absolute coefficient C Di =~^ 
Parasite drag, absolute coefficient D Dp =^ 


C 


Cross-wind force, absolute coefficient 


iu, 

it 

Q 

n 

R 


a 

e 

«a 

7 


Angle of setting of wings (relative to thrust line) 
Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 


Reynolds number, 


VI 


where Z is a linear dimen- 


sion (e.g., for an airfoil of 1 .0 ft chord, 100 mph, 
standard pressure at 15° C, the corresponding 
Reynolds number is 935,400; or for an airfoil 
of 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 

Flight-path angle 
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DISLOCATION THEORY OF THE FATIGUE OF METALS 

By E. S. Machlin 


SUMMARY 

A dislocation theory of fatigue failure for annealed solid 
solutions is presented. On the basis of this theory , an equation 
giving the dependence of the number of cycles for failure on the 
stress, the temperature , the material parameters , and the fre- 
quency is derived for uniformly stressed specimens. The 
equation is in quantitative agreement with the data. Inasmuch 
as one material parameter is indicated to be temperature- 
dependent and its temperature dependence is unknown , it is 
impossible to predict the temperature dependence of the number 
of cycles for failure. A predicted quantitative correlation be- 
tween fatigue and creep was found to exist , which suggests the 
practical possibility of obtaining fatigue data for annealed 
solid solutions and elements from steady-state creep-rate data 
for these materials. As a result of this investigation , a modifica- 
tion of the equation for the steady-state creep rate previously 
developed on the basis of the dislocation theory is suggested. 
Additional data are required to verify completely the dislocation 
theory of fatigue. 

INTRODUCTION 

The failure of materials under an oscillating stress, the 
maximum value of which is lower than that required to cause 
failure in a static tensile test, is termed “fatigue failure.” 
This type of failure is the limiting factor in the design and the 
operation of many rotating parts. In particular, fatigue 
strength has become a limiting factor in the development of 
the gas turbine for aircraft use. Turbine blades, for ex- 
ample, have in many cases failed as a result of fatigue. 

Two methods of eliminating the problem of fatigue failure 
exist: (1) The design of the parts can be changed to minimize 
the effect of oscillating stresses; and (2) the resistance to 
fatigue failure of the materials used can be increased. A 
phase of the second solution, based on obtaining an under- 
standing of the fatigue phenomena, is discussed herein. 

Although many experimental investigations of fatigue 
have been made, few have been concerned with basic con- 
siderations. The work of Gough and his collaborators 
(references 1 and 2), however, is outstanding in the field of 
fundamental fatigue studies. This work has shown that 
the phenomena of plastic deformation (slip) and fatigue are 
closely related. Some theories have been reported (refer- 
ences 3 and 4) that relate plastic deformation to fatigue 
failure. Because these theories are not of a fundamental 
nature, however, they have not yielded basic knowledge 
of the physical properties of materials that affect fatigue. 

With the development of the theory of dislocations, the | 


understanding of the phenomena of plastic deformation has 
greatly advanced (references 5 to 8). Inasmuch as plastic 
deformation and .fatigue have been shown to be related, an 
investigation of fatigue based on the theory of dislocations 
was undertaken during 1947 and is presented herein. In 
this investigation a physical model, based on the dislocation 
theory, was devised to account for crack growth. This 
model was then used as a means of obtaining an equation 
relating the number of cycles for failure to the appropriate 
variables. The equation so obtained was then subjected to 
various experimental checks. This investigation was part 
of a program conducted at the NACA Cleveland laboratory 
to evaluate the physical properties of heat-resisting alloys in 
terms of physical constants that are easily measurable and 
to make possible the synthesis of compositions and structures 
of heat-resisting alloys better than those currently used. 

SYMBOLS 

The following symbols are used in this report: 
a constant 

c proportionality constant relating r a to <r m /2 
d x distance between atoms in slip direction, centimeters 
exp base of Napierian logarithmic system raised to power 
in parentheses following exp 

f fraction, experimental value= 0.374 (reference 8) 

A F g free energy of activation involved in generation of a 
dislocation, ergs per molecule 
h Planck's constant, 6.62 X10 - * 7 erg seconds 

AH g heat of activation per molecule involved in generation 
of a dislocation, ergs 

k Boltzmann's constant, 1.38X10" 16 ergs per molecule 

per °K 

L distance between imperfections in single crystal, of 
order of 1 micron 

M amount of crack growth per crack source necessary for 
failure measured in units of number of interatomic 
spacings 

N number of cycles for failure 

P period of cycle of stress, seconds 
q stress-concentration factor 

q' cq 

H d rate of generation of positive or negative dislocations 
at stress concentration of internal surface in crystal- 
line material, number of dislocations per second 
R g net rate of growth per crack, number of atomic spac- 
ings per second 
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AS g entropy of activation per molecule involved in gen- 
eration of dislocation, ergs per °K 
T absolute temperature, °K 

t time, seconds 

u' shear rate, 1 per second 

Ui antilog of intercept of plot of log u' against a at <r=0 
V atomic volume, cubic centimeters 

x ratio of (h to interplanar spacing of slip planes, (V3/2 
for face-centered and body-centered cubic lattices) 

y mathematical variable 

6 2inxit 

a tensile stress, dynes per square centimeter 
cr m maximum tensile or compressive stress of cycle, 
dynes per square centimeter 

r resolved shear stress, dynes per square centimeter 

t u average resolved shear stress in polycrystalline speci- 

men operating for failure, dynes per square centi- 
meter 

n resolved shear stress due to internal source, dynes per 
square centimeter 

T m resolved maximum shear stress corresponding to a m , 
dynes per square centimeter 

t s average shear stress at sources of dislocations, dynes 
per square centimeter 

a? frequency of cycles of stress application, 1 per second 

CRACK MECHANISM 

Fatigue failure is usually associated with the propagation 
of a crack. This fatigue crack must be formed either as a 
result of internal tensile stresses that exceed the tensile 
strength of the material or from pre-existing submicroscopic 
cracks that grow under the influence of the cyclic stress. 
The assumption will be made that such submicroscopic 
cracks do exist. It is probable that some of the sources of 
stress concentrations in metals are in the form of submicro- 
scopic cracks that occur at the boundaries of mosaic blocks 
(reference 6). 

A theory that would explain the growth of these submicro- 
scopic cracks as a result of the continual reversal of stress 
might provide a basis for the development of a quantitative 
theory of fatigue. The analysis of the mechanism of crack 
formation developed herein is based on the dislocation theory 
discussed in references 5, 6, 8, and 9. A short resume of the 
dislocation theory follows: 

A dislocation consists of a stable arrangement of atoms 
such that, in a region of a few atomic distances, n+1 atoms 
in the slip direction face n atoms across the slip plane. When 
the n+1 atoms are above the n atoms, the configuration is 
called a positive dislocation; when the n -\- 1 atoms are below, 
it is called a negative dislocation. The net result of a posi- 
tive dislocation moving completely through a specimen from 
left to right, or a negative dislocation moving completely 
through the specimen from right to left, is a translation of 
the material above the plane between the n+l atoms and the 
n atoms, with respect to the material below this plane, by 


one atomic distance to the right. If the positive dislocation 
moves from right to left or the negative dislocation moves 
from left to right, the opposite translation takes place. 


Edge of crack after 
partial translation of 
lattice along slip planes 



dislocations 


Figure 1 .— Single crack generating source of dislocations. 

The most probable sources of generation of dislocations 
are stress concentrations. Such stress concentrations occur 
at the ends of cracks. For example, an ellipsoidal crack 
shown in cross section in figure 1, has two points of stress 
concentration (A and B in fig. 1), that exceed the stress 
concentration at any other point along the crack circum- 
ference (reference 10). Dislocations will usually be gen- 
erated at these stress concentrations upon the application 
of a shear stress. 

In the development of a theory for the growth of sub- 
microscopic cracks, an equation that was in agreement with 
the available data was obtained by making certain assump- 
tions about the physical model. Although each of these 
assumptions may have a physical basis, such a basis was 
not immediately evident. The development of the following 
theory therefore presents these assumptions without any 
attempt to justify them physically. 

Cracks similar to the one diagrammatically shown in figure 1 
are assumed to exist in such a manner that points A and B 
act as sources of generation of dislocations. It is further 
assumed that only positive dislocations can be generated 
at point A and only negative dislocations can be generated 
at point B. If the line joining A and B is at a small angle to 
the slip plane, if the generated positive dislocations move to 
the right at point A and if the generated negative dislocations 
move to t he left at point B as a result of the given shear stress, 
the crack will grow. An intermediate position of the right 
side of the crack after the crack has grown to some extent is 
shown by the dashed line in figure 1. 

If all the dislocations that were generated at A and B 
during one half-cycle of stress had returned to their respective 
sources during the next half-cycle of stress, the crack would 
not have grown at the end of a complete cycle of reversed 
stressing. In order for the crack to grow, some of the dis- 
locations generated at the crack sources during the growth 
part of the stress cycle therefore must either disappear from 
the specimen or reach a point from which their return to 
the source is prevented. No dislocations generated at A 
and B in the second half-cycle of stress can enter the crystal 


DISLOCATION THEORY OF THE FATIGUE OF METALS 


3 


lattice because, according to the assumption, only positive 
dislocations are generated at A, only negative dislocations 
are generated at B, and the applied stress during this half- 
cycle of stress would make any dislocations generated at 
these sources return to their respective sources. This 
assumption does not mean, however, that other sources, at 
which none of the dislocations there generated move into the 
crystal lattice during the first half-cycle, will not generate 
dislocations that move into the crystal lattice in the last 
half-cycle. 

If one submicroscopic crack were to account for the failure 
of a macroscopic* specimen, it is probable that more than the 
observed amount of plastic flow would have to occur. If 
failure were to occur as the result of the growth of submicro- 
scopic cracks, it would evidently take place because many 
cracks contribute to the failure. A simple configuration of 
cracks of the type illustrated in figure 1 that can lead to 
fatigue failure with little plastic flow is shown in figure 2(a). 


B 


Operating shear stress 


B 




Operating slip planes 


(a) 


B 



(b) 


(a) Leading to fatigue failure. 

(b) Leading to growth of cracks after partial slip has occurred. 
Figure 2. — Submicroscopic-crack configuration. 


This configuration, after all the cracks have grown to a 
certain extent, would appear as shown in figure 2(b). In 
order that the cracks in this configuration may grow, the 
assumption becomes necessary, as before, that at points A 
only positive dislocations are generated and at points B only 
negative dislocations are generated. The result of positive 
dislocations moving to the left and negative dislocations 
simultaneously moving to the right during a cycle of stress 
favorable for this motion is that relative translations of the 
lattice occur that yield larger cracks than existed prior to 
the translations. Inasmuch as these dislocations disappeared 
from the lattice at points C and D in figure 2(a), the applica- 
tion of a stress opposite to the previous stress would not cause 
any relative translations of the lattice between the regions 
separated by the slip planes A-C and B-D, unless positive 
and negative dislocations were generated at C and D, re- 
spectively. No dislocations are assumed to be generated at 
points C and D. Each crack in the configuration will there- 
fore grow as long as positive dislocations generated at A 


and negative dislocations generated at B disappear from the 
lattice at C and D, respectively. 

In the foregoing discussion, each crack was assumed to 
extend completely through the specimen in the direction 
perpendicular to the plane of the figure. Because the cracks 
were assumed to exist in the mosaic-block boundaries, it 
seems probable, however, that the crack width in the direc- 
tion considered is not more than a few mosaic-block lengths, 
that is, about 10“ 4 centimeters. In order to account for the 
growth of short cracks, it thus becomes necessary to assume 
that such short dislocations may exist and that relative slip 
between mosaic blocks may easily occur. Seitz and Read 
in considering the phenomena of slip with respect to dis- 
locations found that these assumptions would help explain 
the experimental slip phenomena (reference 6, pt. II). 

One of the consequences of the crack-growth mechanism 
just described is that the cracks will begin to grow along slip 
planes. This phenomenon has been observed by a number of 
experimenters (for example, references 2 and 11) on poly- 
crystalline as well as single-crystal specimens. It was noticed 
that failure began to occur along the slip plane and then 
proceeded in the direction normal to the maximum tensile 
stress. 

NUMBER OF CYCLES FOR FAILURE 


The general equation for the rate of generation of disloca- 
tions R d of one type at a stress concentration of an internal 
surface in a crystalline material is given in reference 8. The 


equation, for large values of 


2 qVxfr 
kT ’ 


can be written as 


jy kT 
Rd=-1~ exp I 


—A F\ | 2 qVxff \ 


kT 


kT ) 


( 1 ) 


The net rate of growth per crack R g is the rate of genera- 
tion of dislocations during the growth cycle: 

The stress at a generating source may consist of the applied 
stress and any internal stress that may be present. As a 
result, 

T s =qr=q(r m sin 0+n) (3) 

If it is assumed that the amount of crack growth per crack 
source required for failure M is a constant for all specimens 
and materials and that the stress-concentration values are 
the same for every crack associated with failure and do not 
change with the number of cycles of stressing, the depend- 
ence of the number of cycles for failure on the maximum 
resolved shear stress in the slip plane and direction con- 
sidered is 

p 

M=aT 2 R g dt 


( 4 ) 
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because Z? 8 is independent of N. Substituting for R * in equa- 
tion (4) yields 


M=N^[exp (n^)]J 0 2 «3> ( 2(L It T ) 


(it 


(5) 


For annealed^metals^the'internal' stress r t is* approximately 
zero. Thus, 

T=T m Sill 6 

Changing the variable of integration in equation (5) yields 


, f NkT F / — AFA“1 j ,7r ( 2 qVxf r m sin 


(/« (0) 


The expression log e ( exp (ay sin 0) can be approxi- 
Jo 

mated for large values of ay, by the product ay, as shown 
in appendix A. When this approximation is used in equation 
(6) and equation (6) is solved for log N 


log N i 


dog 


/2t uhM\ A F g 2qVxfT m 

\kT A 2.3 kT 2.3kT 


(7) 


In order to apply equation (7) to polycrystalline speci- 
mens, some assumptions must be made concerning the 
method of fatigue failure of polycrystals. The orientations 
of the crystals comprising the fatigue-failed regions of a 
polycrystal are probably random. The maximum resolved 
maximum shear stress, if each of these crystals was con- 
sidered apart from the others, would therefore vary from 
some minimum (not zero for face-centered cubic or body- 
centered cubic lattices) to the maximum aj 2. Inasmuch 
as these crystals are not separated by stress-free regions, 
internal stresses would be set up between the differently 
oriented crystals. As a result, steady-state slip would occur 
at some stress intermediate to the minimum and maximum 
value of the maximum resolved maximum shear stress. The 
value of this average resolved shear stress should then be 
substituted in equation (7) for polycrystalline specimens. 
If the average resolved shear stress is related to the maximum 

resolved shear stress by the equation r a =c ™ and if q ' =cq, 


then, 


log iV=log ( 


/ 2iro)hM\ 

kT~ ) 


+ 


al 

2.3 kT 


q I xj <j m 

~2wr 


(8) 


Equation (8) applies only for stresses sufficiently above the 
endurance limit (where the log N varies linearly with <r m ) 
because the conditions that determine the endurance limit 
have not been considered in the development of the theory. 
Also, because the assumption was made that a m is constant 
over the cross section of the specimen, equation (8) will 
apply to axial-stress-type loading only. The cases of bending- 
type stresses and torsional stresses are theoretically 
capable of being treated in a manner similar to the case of 
uniform axial stress; however, no attempt has been made to 
evaluate equations of the same type as equation (8) for these 


According to equation (8), the intercept of the straight line 
obtained by plotting log N against a m at a m =0 is equal to 



2i rwIiAI 
kT 



2.3 kT 


If all the parameters in this expression are known, it should be 
possible to obtain a calculated value of the intercept to com- 
pare with the experimental value. The parameters to and T 
are known from the experimental conditions. The param- 
eter M is of the order of magnitude of a mosaic block 
length, measured in units of interatomic spacings, that is, 
about 10 4 . The parameters k and h are known constants. 
The only factor that remains to be determined is A F g . 
A method of calculating A F g by using creep data is pre- 
sented in reference 8. As shown in appendix B, however, 
A F g cannot be completely determined from creep data 
because an assumption made in reference 8 is incorrect. The 
analysis presented in appendix B does indicate, however, 
that a set of fatigue and creep data for the same material 
can be used to calculate a value of the unknown factor. If 
the value of this factor is assumed to be independent of the 
material variable, it can be used to help calculate values of 
A F g using creep data. A table of such calculated values 
using data for copper (fig. 3(a)) to evaluate the unknown 
factor is presented in appendix B. 



(a) Copper. 

(b) Armco iron. 

Figure 3. — Variation of stress with number of cycles for failure at room temperature. 
Frequency, 2200 cycles per minute. (Data taken from reference 12.) 


cases. 
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A check of equation (8) can be made by comparing the 
experimental value of the intercept of the straight line 
obtained by plotting log A 7 ' against <r fri with the value obtained 
from calculation using the term derived from equation (8). 
Data for Armco iron, obtained from reference 12 and plotted 
in figure 3(b), can be used to make the suggested comparison 
of experimental and calculated values of the intercept being 
considered. The experimental value of the intercept deter- 
mined by extrapolating the best straight line through the 
points of figure 3(b) is equal to 22. * The approximate 
uncertainty in this value is ±2. The calculated value is 


log 


"2™AM\ A F g ( 2wX 6o°X 6 -62X10-- ,7 X10'\ 

^ kT / + 2.3 kT og \ 1.38X10 -l6 X300 / 


+ 


2.92X10- 12 


2.3 XI -38 X 10~ 16 X300 


= 24.2 


where o>=2200 cycles per minute and T^SOO 0 K, as given by 
reference 12. The agreement between the predicted value 
of 24.2 and the observed value of 22 ±2 tends to verify 
equation (8). 

An additional check of the validity of equation (8) could 
be made if the value of the stress-concentration factor q 
were known. Because the value of q is unknown, the experi- 
mental value of the slope of the linear portion of the plot of 
log N against a m can be used to obtain a value of q'. This 
calculation is 


2 ' 


2.3 iTX (experimental slope, sq cm/dynes) 
VXxXf 


_2.3X1. 38X300X10 16 ( 0 . 755 x 2240 X 6 . 9X 10 4 ) 

1 1.7 X10 -24 X 1-23X0.374 

= 152 

A similar procedure was followed to obtain a value of q f 
for annealed copper using the line shown in figure 3(a) to 
define the slope. A value of ^' = 164 was obtained. 

The fatigue theory can be further experimentally checked. 
If equation (8) is solved for the parameter Va m , the following 
equation is obtained: 


where 


V<r m =U+WAF g 


(9) 


W= 


a'xf 


Data to check equation (9) were obtained from a number 
of sources. The details concerning the sources and precision 
of data for AF g are contained in appendix B. 

Because the data for the value of a m that corresponds to 
failure in 3.3 X10 6 cycles were difficult to obtain accurately, 
the following table lists, wherever possible, a number of 
sources for a given material: 


Material 

<r m , dynes/sq cm 

Reference 

Aluminum 

3.5X108 

13 

Silver" 

6.2 

14 

Copper 

8.2 

12 


8.65 

15 

Armco iron 

1.74X10 9 

12 


1.90 

16 

Nickel 

2.06 

17 


2.98 

15 


“Data obtained from experiment applying torsional .stresses to a single crystal. Because 
the corrections involved are in opposite directions, the value given is a fair approximation. 


Inasmuch as different sources yielded different values of 
the stress corresponding to a certain number of cycles for 
failure for a given material and because the values of AF g 
were obtained by extrapolation, the data were plotted as 
regions of approximately equal probability corresponding to 
each material. The data are therefore plotted as rectangles 
or lines in figure 4. The data in this figure correspond to 



Figure 4.— Dependence of stress corresponding to failure in 3.3X10* cycles on material 
parameters at room temperature. 


failure in approximately 3.3 X10 6 cycles; the frequency 
used to calculate U is 2200 cycles per minute; the tempera- 
ture is 300° K. 

The theoretical curve corresponding to these data is, ac- 
cording to equation (9), a straight line of slope W and inter- 
cept U. The value of q ' was considered to be 164, as pre- 
viously found for annealed copper. Because all the other 
factors involved in the calculation of U and W are known, 
calculation of the theoretical values of U and W is possible. 
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The calculated values of V and W are — 157X10" 16 ergs 
per molecule and 0.0132, respectively. The theoretical curve 
corresponding to these values of U and W is plotted as the 
straight line in figure 4; the theoretical line is in agreement 
with the data. Figure 4 represents a correlation between 
creep and fatigue, inasmuch as the ordinate is a parameter 
determined from fatigue data and the abscissa is a parameter 
determined from creep data. 

Two other independent checks of the theory are possible: 
The first is a check of the temperature dependence of the 
number of cycles for failure; the second is based on the de- 
pendence of the number of cycles for failure on frequency. 

The assumption is made that q is independent of temper- 
ature. If equation (8) is solved for a m , the following equation 
is obtained: 

* m =A+B(C-log T)T ( 10 ) 

where 



Data for Armco iron obtained from reference 16 are shown 
in figure 5. A cross plot of these data (temperatures in °K) 
is shown in figure 6, which also gives the theoretical curve 
corresponding to this experimental plot. Data used to 
compute the theoretical values were obtained as follows: 
Values of A H g /V correspond to values of A\V from figure 6 
of reference 8; values of A S & were taken from the values 
listed in the table given on page 16 of reference 8, corrected 
to account for the factor determined in appendix B; the 
value of q' was taken as 164; the frequency w was taken as 
40 cycles per second, as reported in reference 16. 


.c 



Figure 5. — Temperature dependence of curve for c m plotted against N for Armco iron. 
(Data taken from reference 16.) 


A study of figure 6 indicates that t lie theoretical curve 
yields too small a temperature variation of a m . Therefore, 
equation (10) apparently does not agree with the experimental 
data. Two explanations of these results are possible. The 
first explanation is that the data may be in error. The large 
scatter for the data corresponding to 500° C, as shown in 
figure 5, indicates that this explanation may be true. Fur- 
thermore, meaningful temperature studies on direct-stress 
specimens are difficult to perform because the specimen may 
easily become misalined, owing to the use of long specimen 
holders. The misalinement may then be a function of 
temperature and yield an apparent temperature effect, as 
shown in figure 6. 

The second and more probable explanation is that the 
stress concentrations associated with the cracks in the mosaic- 
block boundaries are temperature-dependent and, hence, the 
assumption made in plotting the theoretical curve in figure 
6 is not valid. In this case, predicting the temperature 
dependence of a m becomes difficult inasmuch as no apparent 
basis exists for determining the dependence of the stress- 
concentration factor on temperature. 

A final check of the dislocation theory of fatigue can be 
obtained by determining whether the dependence of the 
number of cycles for failure on the frequency, as given by 
equation (8), is experimentally obtained. Equation (8) also 
indicates that log N would be linearly related to log to 
at constant maximum stress and temperature. The experi- 
ments reported in reference 18 yield results that qualitatively 



Temperature , °K 

Figure 6.— Temperature dependence of stress corresponding to failure in 10 s cycles for 

Armco iron. 
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verify the prediction. In reference 18, direct-stress experi- 
ments performed at 7000 cycles per minute required a 
greater number of cycles for failure than tests at 1200 cycles 
per minute. The tests at 1200 cycles per minute, however, 
were performed on a different type machine than the tests at 
7000 cycles per minute; therefore, the results would not be 
expected to agree quantitatively. The increase in log N was 
nevertheless sufficiently large that it is improbable that the 
difference in machines could be solely responsible for the 
effect noticed. Reference 19 indicates that no frequency 
effect occurs at low values of frequency. An examination of 
the data, however, reveals that more data would be required 
in order to be able to draw any conclusion concerning this 
effect. 

EVALUATION OF RESULTS 

Sufficient data to verify equation (8) completely are un- 
available. Inasmuch as this equation is in quantitative 
agreement with available data, further experimentation 
designed to verify the dislocation theory of fatigue is justified. 

The theory was developed for annealed elements or solid 
solutions; the data used to check the theory were obtained 
only from annealed elements. It seems reasonable, however, 
that data for annealed solid solutions would also be in 
agreement with the theory. 

It is improbable that equation (8) would apply to the case 
of precipitation-hardened or strain-hardened materials be- 
cause certain modifications, such as the introduction of an 
internal-stress term, would be necessary. This subject is 
sufficiently complex to justify an independent investigation. 


The uncertainty in the data appears sufficiently large to 
allow reasonable doubt in the constancy of the stress- 
concentration factor q for all annealed elements and solid so- 
lutions. The possibility exists that an experimental plot of 
q'Va m against AF g would also be linear; however, the data 
are insufficient to resolve this question. 

The number of cycles for failure at some given stress and 
temperature have been experimentally observed to follow a 
statistical distribution (reference 20). An explanation of this 
phenomenon on the basis of equation (8) would probably 
require that one or more of the following factors, M, AF g) and 
q assume a range of values for a given material. It seems 
reasonable that the factors M and q may vary from specimen 
to specimen inasmuch as these quantities are probably de- 
pendent on manufacturing variables such as casting, drawing, 
and heat-treating. On the basis of the dislocation theory, 
accounting for the scatter of data normally obtained in 
fatigue tests is thus possible. 

In figure 4, which helps substantiate equation (8), a corre- 
lation between fatigue and creep is established as predicted 
by the theory and fatigue failure is proved to depend on 
plastic deformation. With the assumption of a constant value 
of q at room temperature, the prediction of fatigue data from 
a knowledge of creep data for annealed elements and solid 
solutions is made possible for practical purposes. 


National Advisory Committee for Aeronautics, 
Flight Propulsion Research Laboratory, 
Cleveland, Ohio, September 12 , 1947. 


APPENDIX A 

EVALUATION OF DEFINITE INTEGRAL 


The definite integral 

exp {ay sin 6) dO 

was evaluated for a series of values of ay by numerical in- 
tegration. The value of a was set equal to 10.56X10" 8 to 
correspond approximately to the expected values of the 
physical quantity it represented. A plot of the log of the 
integral against y is shown in figure 7. For values of y >10 8 , 
the equation of the curve is 

logj^ exp (10.56X10 8 2/ sin 6) dd^by 

where b is the slope of the line in figure 7, 4.44 X10~ 8 . As 
a result, 

loge J exp (10.56X 10~ s 2/ sin 6) d0^l0.22yX10~ 8 


or 



exp {ay sin 6) dd^ay 
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APPENDIX B 

EVALUATION OF AF t 


The equation 


tt '= 10- 4 


kT 

T exp 


( ~ AF t \ 2 g Vxfr \ 
\ kT "*■ kT ) 


obtained from reference 8 was first used to calculate values of 
AF g from steady-state creep-rate data. This calculation was 
accomplished by plotting log u' against a, setting the experi- 
mental value of the intercept at <r =0 equal to 



kT AF g \ 
h 2.3 kT) 


and solving for A F e because all the other parameters are 
known in the expression for the intercept. When these 
values of AF S were used to calculate the intercepts of the 
curves of log N plotted against <r m and Va m plotted against 
AF e , the calculated intercepts were found to exceed the 
experimental intercepts by a factor that varied between 
1X10 ~ 13 to 3X10 -13 ergs per molecule in the value of AF g . 
In attempting to find an explanation for this discrepancy, 
an assumption made in reference 8 regarding the value and 
constancy of the quantity 10~ 4 in the preceding equation 
was found to be erroneous. The value of 10“ 4 was used as 
an approximation for d^L. In the derivation of the fore- 
going equation, the value of lO -4 was based on an assumption 
that one source of generation of dislocations exists per 
mosaic block. Thus, according to this assumption, the 
strain due to the passage of one dislocation through a mosaic 
d 3X1 0 -8 

block is 3 y jq^ 4 ~ 10 -4 . Evidence exists, however, for 

single crystals, which indicates that the spacings between 
slip planes may be many multiples greater than 3X10 -4 . 
For example, figure 8 , which was taken from reference 21 , 
indicates that the spacing may be as large as 0.3 centimeter 
find figure 9 indicates that 1/L (number of slip planes/cm) 
may vary between 0 and 7500. Thus, the assumption that 
L is a constant for polycrystals of magnitude 3X10 -4 is 
probably incorrect. It remains to be shown, however, how 
L may be calculated. 

Tempera ture 
(°C) 



Figure 8.— Pure aluminum sample elongated at different temperatures by application of 
stress. Most of plastic flow occurs along slip bands. Not etched. XI 5. (Taken from 
fig. 14 of reference 21.) 


Without a physical model of the mechanism of strain in a 
polycrystal, it seems that the simplest assumption is to leave 
the factor dJL as an unknown constant 10 c . The factor C 
can then be determined from a set of fatigue and creep data 
for the same material as follows: From equation ( 8 ), it is 
evident that the intercept of the straight line obtained by 
plotting log N against <r m at <r w = 0 is 

, /2wuhM\ . AF g 

'°z{-Tr-)+2WT 

The data for annealed copper shown in figure 3 (a) were used 
in this calculation. The data for dry purified air were used 
to locate the slope of the line taken through the points for 
the determinations in vacuum. The intercept of this line 



Figure 9.— Stress dependence of number of slip planes per millimeter. (Taken from fig. 24 

of reference 26.) 
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was found to equal approximately 15. Hence, 


&Fg— 


( 22 
2irX ~6( 
L 


2200 

60 


X6.62X 10 _27 X 10 4 


38X 10 -16 X300 


2.3X 


1.38X10" 16 X300 


= 2.01X10~ 12 ergs/molecule 

where co=2200 cycles per minute and 7=300° K, as given 
by reference 12. 


From the creep equation, 


log u' = C+log 


kT A F g . 2 qVxfr 

h-2WT + kT~ 


it is apparent that the intercept log u t of the straight line 
obtained by plotting log u' against rat r=0 is 

, ^ ! kT A F g 

log Ui=C+ log T 

An estimate of log u t at 300° K can be obtained from data 
given at other temperatures by plotting (t log against 

yp> because A F g is linearly dependent on T , and by extrapo- 
lating the best straight line thus obtained to 7=300° K. 
From the data of references 22 and 23, an average value of 
— 15 was obtained for log u t . When this value was substi- 
tuted for log u t and A7*=2.01 X 10“ 12 in the preceding equa- 
tion, solving for C yielded 

C=- 6.7 


For the purpose of obtaining values of AF g from creep data 
for other materials, C was assumed to be a constant. When 
this value of C was used, values of A F g were obtained from 
creep data and are given in the following table together with 
a list of the sources of the respective data: 


Material 

A F 0 , ergs/molecule 

Reference 

Aluminum 

1. 57-1. 70XKH 2 

21,24 

Silver - 

1.85-2. 26 

21 

Copper 

1.94-2. 12 

22,23 

Armco iron - 

2. 74-3. 09 

25 

Nickel 

3. 13-3. 63 

21 
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Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Designation 

Sym- 

bol 

Longitudinal 

X 

Lateral 

Y 

Normal 

z 




Force 
(parallel 
to axis) 
symbol 


X 

Y 

Z 


Moment about axis 

Angle 

Velocities 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 

Rolling 

L 

Y >Z 

Roll 

d> 

u 

X) 

Pitching 

M 

Z >X 

Pitch. 

o 

v 

r 

0 

Yawing 

N 

X >Y 

Yaw 

* 

w 

r 


Absolute coefficients of moment 

C — L C — ^ 

( rolling) (pitching) 


o-JL 

Cn ~qbS 

(yawing) 


Angle of set of control surface (relative to neutral 
position), 8. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 


D 

Diameter 


p 

V 

Geometric pitch 

P 

Power, absolute coefficient Cp= — ^ 
1 pn z D b 

P/D 

Pitch ratio 

c, 

Speed-power coefficient 

V' 

Inflow velocity 

V, 

Slipstream velocity 

rn 

V 

Efficiency 

T 

Thrust, absolute coefficient C T — — 

pn*D* 

n 

Revolutions per second, rps 

Q 

Torque, absolute coefficient C 0 ~ 

pn z Lr 

$ 

Effective helix angle=tan _1 ( ) 


5. NUMERICAL RELATIONS 


1 hp=76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower=0.9863 hp 
1 mph=0.4470 mps 
1 mps=2.2369 mph 


1 lb=0.4536 kg 
1 kg=2.2Q46 lb 
1 mi= 1,609.35 m=5,280 ft 
1 m=3.2808 ft 


